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We study the lepton-fiavor-violating reaction /i + e~ — > e + e _ within two extensions of the standard 
model that include heavy neutrinos. The reaction is studied in the low energy limit in the form of 
muonium decay M — > e + e~ and in the high energy regime of a muon collider. The two theoretical 
models we consider are: model I, a typical see-saw model that violates lepton flavor and number by 
inclusion of extra right handed neutrinos, and model II, a variant where lepton number is conserved 
and which includes extra right handed as well as left handed neutrinos, singlets under the gauge 
group. We find for muonium decay into e + e~ the extremely small result Br(M — > e + e~) < 10 -19 in 
both scenarios. Alternatively, for [i + e~ collisions up to y^i ~ 50 GeV we find cr(/i + e~ — > e + e~) < 
1(T 5 fb, while for energies above the W + W threshold we find a(fi + e — » W + W ) up to ~ 1 fb. 
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I. INTRODUCTION 



A current issue in particle physics is the understanding of the neutrino sector. It is by now consistent with 
experimental results to assert that at least some of the standard neutrinos must have small (but non-zero) masses 
and these mass states should mix under the weak interactions As a consequence, the lightness of the standard 
neutrinos {jn v < 1 eV) compared to the charged leptons (> MeV) remains to be understood. If neutrinos are of a Dirac 
nature, non-zero masses can appear in the standard model (SM) by inclusion of (sterile) right-handed neutrinos [2J. If 
neutrinos are of a Majorana nature -a more general case-, more appealing solutions exist in the context of extended 
gauge structures. An interesting solution is provided by the so-called seesaw mechanism within SO[VS) or left-right 
symmetric models of interactions. In conventional seesaw models, the effective light neutrino masses are within the 
scales of eV to MeV via a relation involving the hierarchy between very large Majorana masses and Dirac masses 
comparable to those of the charged leptons Q . Another possible solution was investigated in the framework of heterotic 
superstring models Q with Eq symmetry or certain scenarios of 50(10) models [5], where the low-energy effective 
theories include new left-handed and right- handed neutral isosinglets and assume conservation of total lepton number 
in the Yukawa sector. 

One avenue to study the neutrino sector is to consider lepton flavor violating processes (LFV) of standard particles, 
which are absent in the SM and depend on neutrino masses and mixings. Free lepton decays like fj, — > ej n — > 3e, 
r — > 3/i, etc. are among the most direct LFV processes and they have been studied in detail @-[IH,[l^]- Other ways 
like muon-to-electron conversion in nuclei and muonium-to-antimuonium conversion are also possibilities [l2l . [T3 . [lq | . 
Yet another possibility is to consider lepton- lepton scattering or bound state decays (e.g. muonium decay). These 
modes contain the same underlying physics as free lepton decays, but differ in the experimental conditions and, 
from the phenomenological viewpoint they are sensitive to other combination of parameters because of their different 
kinematical regime. 
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In this work we study the process /i + e~ — > e + e~ at low energy in the form of muonium decay (muonium is a 
non-relativistic /i + e~ bound state) and at high energy in the form of fi + e~ collisions. Muonium is formed when a 
/i + slows down inside a material and captures an electron, forming a bound state due to Coulomb attraction. The 
fraction of muons that actually end up as muonium depends very much on the material, ranging from small fractions 
in some to nearly 100% in others. Most muonia decay as M — > e + e~ u^Ue, which is simply a muon decay with an 
electron as spectator [lij]. A more rare mode, albeit standard, is annihilation M — ► v^Ve, whose rate is ~ 10~ n times 
smaller. Here we study the mode M — > e + e~, which violates lepton flavor and is thus forbidden in the SM. From 
our point of view, muonium decay is interesting because it is a non-relativistic [i + e~ process, and thus depends on 
different combinations of parameters than in high energy processes. Complementary to it, we also study the collision 
(i + e~ — > e + e~ at high energy. This process clearly requires a muon collider. We will not be exhaustive in the study 
of new physics at muon colliders but focus on this particular channel within the LFV models previously mentioned. 
In section II we give a brief review of the seesaw type models in question; in section III we describe the generic 
amplitudes for /i + e~ — > e + e~ processes within these models; section IV contains the application to muonium decay, 
section V the application to high energy collisions and section VI summarizes the results and conclusions. 



II. TWO MODELS FOR HEAVY NEUTRINOS 

Here we will briefly summarize the two models under consideration. Further details were given in previous works (TH 
[l6j . The main features of these models are (i) they contain heavy neutrinos, (ii) in the mass eigenstate basis, the heavy 
neutrinos couple with charged leptons of mixed flavor under weak interactions (just like Cabbibo-Kobayashi-Maskawa 
mixing in the quark sector), and (iii) the neutrinos are generally of Majorana type. 

Model I: This is a usual seesaw model [3|: it is the standard model with Nl generations (i.e. left-handed neutrinos 
UiL) and an equal number Nh = Nl of right-handed neutrinos vm, singlets under the gauge group. The neutrino 
mass terms, after gauge symmetry breaking, are: 



i(^W~0 + h.c. (D 



- C v 

mas„ ^ ^ it , j \ VR I 

—T 

The superscript c at some of the fields denotes charge-conjugated fields ip c = Cip , where C = —i-f 2 ^ in the Dirac 
representation. The 2Nl x 2Vi-dimensional matrix M. has a seesaw block form: 

. / m D \ fab \ ( M x \ . . 

M =\ m T m „ ), m D =( _ iSa A ), m M =( M „ ) > (2) 



The (real) parameters a, 6, c, d in the Dirac mass matrix mo are smaller than the Majorana mass parameters Mi 
and M2 (M2 > Mi > 100 GeV) by at least an order of magnitude. 

The squared ratio between the Dirac mass (mo) and the Majorana mass (wim) scales gives the order of magnitude 
of the heavy-to-light neutrino mixings (s^ 1 ) 2 = J^h \Um\ 2 ~ |w£>| 2 /|mM| 2 - The physical light neutrino masses are of 
the order ofm l/[ight ~m D /mM- The very low experimental bounds m Vl . M ~ 1 eV impose in Model I severe constraints 
on the \mjj\ -C |mjv/| hierarchy required. The present experimental bounds on the heavy-to-light mixing parameters 
( s l) 2 ~ I to d| 2 /|^m| 2 (~ 10~ 2 , see below) present another set of constraints on the model. 
Model II: This model contains an equal number Nl of left-handed (Sn) and right-handed (^ir) neutral singlets 
Furthermore, the form of the mass matrix M. leads to the conservation of total Lepton Number, although lepton flavor 
mixing is still possible. The neutrino mass terms, after electroweak symmetry breaking, have the form 
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, , , m D 

v R \ Sl) M\v r \+ h.c. , M = \m T D m T M , (3) 



1 /_ ^- . . I ~ L 

VL 



/ V m M 



and ijid and tum are given in Eq. ©. The mass matrix Ai is (Nl + Nh) X (Nl + V^)-dimensional, Ng = 2Nl 
(the Dirac block is Nl X Nl -dimensional) . The model, in its unperturbed form, predicts for each of the Nl 
generations a massless Weyl neutrino and two degenerate Majorana neutrinos [ijidli (^Vff = 2iVx). Therefore, the 
seesaw-type restriction m Ulight ~ m 2 D /mM of Model I is not present in Model II in its unperturbed form. However, 

present experimental bounds on heavy-to-light mixing parameters ( s L 1 ) 2 ~ \m D \ 2 /\m M \ 2 (~ 10" 2 ) do impose a 
certain level of hierarchy \mr>\ < between the Dirac and Majorana mass sector, but is in general significantly 

weaker than in Model I. Nonzero masses of the Nl light neutrinos can be generated in Model II by introducing small 
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perturbations in the lower right block of M, i.e., small Majorana mass terms for the neutral singlets SiL, and this 
does not significantly affect the mixings of heavy-to-light neutrinos. 

In either model, the mass matrix M can always be diagonalized by means of a congruent transformation involving 
a unitary matrix U 

UMU T A* = M d , (4) 

where Md = diag(mi, 7712, ■ ■ ■) is the nonnegative diagonal mass matrix, and A* is an arbitrary diagonal unitary 
matrix: (A*)jj = &yA* and |Aj| = 1. The Nl, + Nh mass eigenstates rii are Majorana neutrinos, and are related to 
the interaction eigenstates v^l and VjR by the matrices U and A* as: 

Nl+Nh / c \ N l +Nh 

" L ] = E u *- n ^ (It = E u *» A * n « > ( 5 ) 

i=l ^ ' o Z=l 

The first Nl eigenstates n$ (i = 1, . . . ,Nl) are the light partners of the standard charged leptons. The other Njj 
eigenstates are heavy. It can be checked from relations (0 that the relation (n^) c = \*nm holds, and thus is 
recognized as the creation phase factor [l||[2(| of the Majorana neutrino Uj. 

In the mass basis, a Nl x [Nl + iV#)-dimensional mixing matrix B for charged current interactions, and a (Nl + 
Nh) x (Nl + A?#)-dimensional matrix C for neutral current interactions can be introduced: 

n l 

B li = U* l , dj =Y,U ia U* a . (6) 

0=1 

Accordingly, the leptonic weak interactions are (cf. (l6l. l2i|) 

£ lnW (x) = \-ijjL) W-(x) E E ^(^P^-^ + h.c. , (7) 

/ v JVi N L +N H 

\ V W / j — ^ 

C nnZ (x) = f- J Z » E ^Wt" - C IjPr) + h - c - ' ( 9 ) 

N L +N H 

Cnncix) = jTj- G° (x) ~S~] rii(x) {i (mj - rrii)lmCij + 7s(mj + mj)Re(7y} rij (x) , (10) 

where G~ and G° are the Goldstone bosons that appear in a general gauge. Here, the conventions of Itzykson and 
Zuber[22| are used; Pl = (1 — Ys)/2, Pr = (1 + 7s)/2; is the i-th negatively charged lepton (ii = e~, Z 2 = M - )> 3™ 
is the SU(2)l coupling constant (g\, = &GpMyy / 'V2) and nj are the Majorana neutrinos (mass eigenstates). 

In all our numerical calculations we will just use two generations (i.e. Nl = 2), since a third generation does not 
make a sizable change in our results, considering the uncertainties involved. 



III. TRANSITION AMPLITUDE FOR fi+e' -> e+e" 

The amplitude for this process is similar to those of Ref. [7fl for lepton decays. Just as in that case, the amplitude 
for /j + e~ — > e + e~ can be calculated neglecting external momenta and external masses inside the loop, provided the 
neutrinos in the loop are heavy. Let us define the kinematics of the process as: 

M + (fci) e-(fca) -> e+(k 3 ) e -(k 4 ). (11) 

The transition amplitude is a sum of photon- and Z- penguin diagrams (Fig. 1) and box diagrams (Fig. 2). There 
are two types of penguin diagrams in each case: those where the photon or Z is in the t-channel and those where it 
is in the s-channel (see Fig. 1). There is a relative minus sign between those diagrams due to Fermi statistics. The 
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contribution to the amplitude from the photon penguins has two form factors, and G^ M 



iT 7 



. o? w sin 2 6 
2M 2 , 



fy(fci) ( P^ ( 7 p - ^)P £ - iG? 4 ^-(m e P R + m M P L ) ) Ue (fc 2 ) 



u e {k A )^ p v e {kz) 



v e (k 3 ) <-> u e (k 2 ), q -> q 



{q = k 1 + k 2 , q' = ki-k 3 



(12) 



where q u = <7 2 /(47r), Pr/£ = (1 ± 7s)/2 are the chiral projectors, and er po - = (i/2)[7 p , r ) a \. The term g p ^ vanishes 
when contracted to the vector current of the e + e~ pair. Similarly, the Z penguins contain one form factor F^: 



i^L^l F e / 
4M 2 ,} z 



v^k 1 ) 1 p PLU e (k 2 )\ x [u e (fc 4 )7p((l-2sin 2 ^)P L -2sin 2 ^P fl ) We (fc 3 ) 



(13) 



and the box diagrams also one form factor, Fq 



AM, 



F' 

2 B 
W 



v^(ki)j p P L u e (k2) x ?i e (fc 4 )7pP L w e (fc 3 



(14) 




e~(k 2 ) 



e (fc 4 ) 



e+(fc 3 ) 




(a) (6) 

FIG. 1: The one-loop penguin diagrams for /i + (fci)e~(fc2) — > e + (fc3)e~(£;4): (a) s-channel penguin and (b) t-channel penguin. 
The blob represents the loop involving IT, G bosons and neutrinos. The diagram (b), i.e. t-channel, corresponds to the 
contributions in the last lines of Eqs. {12} and (|13[) . 
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e (fc 4 ) 
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i w,g ; 



y[Z (fc 4 ) 



(a) 



(b) 



e~(k 2 ) 



n, 



e (fc 4 ) 



« 1 x 




;W,G 




► 1 x 





e+(fc 3 ) 



e-(fca) 



-x- 



/'W,G 



-x- 



e (/c 4 ) 



e+(A; 3 ) 



(c) 



(d) 



FIG. 2: The box diagrams for fi + (ki)e (/C2) — > e + (hi)e {k^). The diagrams (b) and (d) are obtained from (a) and (c), 
respectively, by the exchange e~(fc2) «-» e + (fca) and by the corresponding change in the flow of time. 



The form factors F*», G e f, F e / and F^ ee are 
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G eAl 



N L +N H _^ 

53 B* ei B, a F^xt) « J] KiB^iF^X!) , (15) 

»=1 />A r t 

^ G 7 (x,) « ^ B^B^G^xj) , (16) 

»=1 />A r i, 
N L +N H 

53 % F z ( Xi ) + iJ z (as* , ccj ) + G* G z (as* , Xj ) (17) 

53 B* eI B^F z ( Xl ) + 53 fl^B^[Cjjffz(a!j,a:j)+C7? J (Gz(a:i,a:j)-Gz(a!/,0)-Gz(0,a:j))] 

/>JVi I,J>N L 

+2 53 B^B^GzfaQ) , (18) 

53 [2|flei| 2 B^ Fb«(x ( , Xj) - (B^B^B^X* G Box (x u Xj )] (19) 
2 53 iB^pB^^j [F Box (xi,xj) + F Box (0,0) - F Bo x(z/,0) - F Box (0, arj)] 

I,J>Nl 

53 (B* eI ) 2 B eJ B^\j\* I G Box (x I ,x J ) + 2 53 B^B |i j(F B «(0,a:j)-FBax(0,0)) , (20) 

I,J>N L J>N L 



where = m|/M^ (ro; is the mass of the Majorana neutrino of type i inside the loop), and the functions Fj(x), 
G 7 (x), Fz{x), Hz(x, y), Gz{x, y), Fb x{x, y) and G Bo x(£, y) appear from the loop integrals 0, in the approximation 
where external masses and momenta are neglected. The explicit expressions of these functions are given in Ref. [lg]. 
Fbox corresponds to the diagrams of Figs. HJa) and (b), and Gbox to those of Figs. [He) and (d) which exist only 
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for Majorana neutrinos. Eqs. (|15[) - (|19p involve summation over all 2-component neutrinos, where Nl is the number 
of standard (light) ones and Nh the number of extra (heavy) ones (= Nl,2Nl in Models I, II, respectively). All 
these expressions are identical to those obtained in Ref. p} for neutrinoless lepton d ecay s, except for the sign in front 
of the term proportional to Gbox(^z,^) _ for details on the latter point, cf. Ref. [16|. Furthermore, in Eqs. (fT5|) - 
(|20P the expressions in a modified form are included which involve only the light-to-heavy mixing elements Bn 
(£ < Nl; I > Nl); they are obtained by taking into account the unitarity of the U matrix, and the fact that the light 
neutrino masses are me <C Myy- 

The parameters Xi are the creation phase factors [20| of the Majorana neutrinos rii and rij , respectively (jij = X*rij), 
which appear in the matrix A of Eq. |QJ . In the case of no CP violation and real mass matrix M. , we can choose the 
unitary matrix U as a real orthogonal matrix which diagonalizes the symmetric and real matrix M. ; then some of the 
Aj's are +1 and others are —1 [16|. Nonetheless we can choose all Xj = +1, but then in general the unitary matrix 
U will not be real. 

The total amplitude for /i + e~ — * e + e~ is then the sum of (fT2)) . (fT3|) and ([M]) . 



IV. THE DECAY OF MUONIUM INTO e~e+ 

The calculation of the decay width of muonium M into e~e + is determined by the LFV transition fi + e~ — ► e + e~ 
at low momentum and by the muonium bound state wavefunction. The formula for the decay width of muonium M 
into a final state / is analogous to the expression for positronium decay: 



dT(M -» /) = 



4fiM 



|(/|2>,e)| 2 dLips. 



(21) 



Here [i and M are the reduced mass and total mass of the bound state, respectively; "0(0) is the spatial wavefunction 
of the bound state at x = 0; the bar over the matrix element means average over the initial spin components; and 
dLips is the usual Lorentz invariant phase space of the final state, dLips = 5 A (pj — pi) d 3 p% d A p±/ {IQt^E^E^). 

The width can also be expressed in terms of the cross section of the constituents at low momentum, er(/i + e~ — > /), 
since the latter involves the same transition element and phase space integral: 



r(M^/) = |v(o)| 2 



'12 



2/iM 



a(u.+e- - /). 



(22) 



Here S12 = 2 Ecu Pcm is the flux factor. At low momenta, the factor 5x2/2/iM reduces to the relative velocity 
between the constituents, u re i. Since muonium is a rather nonrelativistic Coulomb bound state, the wavefunction is 
calculable: 



| V ,( )| 2 = ^L^! 



7T (1+me/mJ 3 



(23) 



and is independent of the final state / into which the constituents disintegrate. On the other hand, the disintegration 
process itself is totally contained in the cross section cr(/u + e _ — > /). Using the two previous expressions, the decay 



M 



e e n 



results in: 



7r (l + m e /m M ) 3 



er(e — > e e + ) 



(24) 



Therefore, the calculation of this decay rate reduces to the calculation of the scattering cross section <r(e /i + 



) 

at low (nonrelativistic) momentum. We remark that the factor v IC \ in Eq. (|2"4")) will cancel the corresponding factor 
1 /v ie \ contained in the cross section. Now, let us find the decay rate within Models I and II described in the previous 
section. 

Although our expressions for the amplitude and form factors are the same as those of Ref. Q for lepton decays, the 
matrix element for muonium decay becomes a different combination of form factors, because the kinematic regime in 
muonium decay is quite different than that in free lepton decay: in muonium, we should take the initial 3-momcnta 
as non-relativistic (or vanishing, since the amplitudes are smooth functions of momenta). The result for the total 



7 



\% 



total | 



amplitude squared, summed over final spins and averaged over initial spins, is then: 

(m e m\ + 2 m e 2 ml + m e 3 m M ) \2Ff + F^ ee | 2 

4 sin 2 M (2m e ml + 3m e 2 m 2 + 3m e 3 m M ) Re [(2F|" + F^ ee )(F^ - F e /)*\ 
12 sin 2 9 W (m e to 3 + 2 m e 2 m 2 + m e 3 m M ) Re [{2F e / + F^ ee )G e ^ *1 
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4sin 4 6>„, (7m e m M 3 + 12m e 2 m 2 + 9TO e 3 rrv) - F^| 2 

4 sin 4 0„, (-2 to 4 + 12 m e m 3 + 36 m e 2 m 2 + 18 m e 3 m M ) Re [(F^ - Ff)G^ 



4 sin 4 9„ 



— + 2 m 4 + 8 m e m 3 + 24 m e 2 m 2 u + 9 m e 3 m tJ 



\ G T\ 



(25) 



From this expression and Eq. (|24[) one obtains the muonium decay rate, after integration over phase space. Since 
the amplitude is independent of the kinematics in the non-relativistic limit, the integration is simply J dLips — 
(l/87r)-v/l — 4m 2 / (m M + m e ) 2 , and we have 

T{M^e+e-) = — '- T1 / 1-4 \'T, 



32n 2 m M (m^ + m e ) 3 



total 



327r 2 m,. 



O 



,2 \ "I 



it; 



total | 



(26) 
(27) 



where |7t t a /| 2 is given in Eq. (f25l) . Now, comparing Eq. (j26|) with the tree-level expression for the muon width 
r(/i — > all)(°' = G 2 7 m^/(1927r 3 ), which is comparable to the muonium width, we obtain the branching ratio: 



Br(M -> e+e~) 



T(M -> e+e~) 

r(^->oH)(°) 



(m M + m e ) 3 



(m M + m e ) 2 



-5 



sin 



red. 



1.77-10" 



B T 



ed. 



where we defined a reduced branching ratio -Bred. 



B re d. = 



4 M^m 



17; 



total 



(28) 



(29) 



Hi' 



with \Ttotai\ 2 given by (|25|) . Provided the neutrino masses are below a few TeV, a rather good approximation for -B re d. 



is obtained by expanding Eq. 

Bred. = 



sin 

1 

4 



to second order in (m e /m M ): 

sin 4 9 w Re [(F 7 - F Z )G* 1 - |G 7 | 2 ] 



9 W | G 7 1 — 2 



771 e 



|F B + 2F Z + 4 sin 2 lu (F 7 - F z ) + 6 sin 2 0„,G O 



4sin 4 6», i , (3|F 7 



|G 7 | 2 



3|F z | 2 -6RcOFzF 7 *)) 



(30) 



For simplicity of notation we have omitted the superscripts e/i and e/iee in the form factors. We see that the box 
diagrams (i.e., the form factors -Fbox) have in general only a small effect on this process. This contrasts with muonium 
to antimuonium conversion, where the box diagrams are the only contribution [16j. This is a bit disappointing, since 
it is only in the box diagrams where the true Majorana character of neutrinos appear, as seen in Figs. 2(c) and (d). 
The rest of the diagrams are sensitive to the neutrino masses and mixings, irrespective of their Dirac or Majorana 
character. In the next section we will see, however, that the free /j + e~ scattering at high energy depends on a 
different combination of form factors, where the box diagrams do enter. In this sense, /i + e~ scattering at high energy 
is complementary to muonium decay, when testing this LFV process. 

In order to obtain numerical results, the parameter space of the models must first be reasonably narrowed down. 
First the parameters appearing in the Dirac matrix mn, Eq. ((2]), are chosen in such a way as to fulfill, at fixed 
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FIG. 3: The branching ratio T(M — > e+e~)/T(fi -* all) 1 - ), Eq. {28JI in Model I, for the (approximately) optimal choice [cf. Eq. ||3TJ] for 
the Dirac mass parameters, as a function of the neutrino mass ratio M2/M1 (> 1), and for different Mi = 100 GeV, 300 GeV, 1 TeV, 10 
TeV, 20 TeV and 30 TeV, as indicated. The last two curves are shown only up to M2/M1 5 and 3, respectively, due to the perturbative 
unitarity bound [7). 



Majorana masses Mi and M 2 , all the experimental requirements, and at the same time to aj 



amplitudes (e.g. |G 7 | etc.) appearing in the decay width of Eqs. 



and 



(see Refs. 



roach maximal squared 
la l for details), in order 



to find the most favorable case to observe such rare processes. This optimal condition fixes the model parameters as 
follows: 



Model I : a 
b 



M lSL 



-d = ay/M 2 /Mi , 5i= S 2 
M 2 s v £ 



tt/2 



Model II : a 

b = c=ax(s u jf;/sl e ) , d = a x (s^'/s^) 2 , S t = S 2 = 



(31) 



(32) 



Here. 



0#) 2 



common value s v £ 



Eh \Bihf (N L <h< N L 
s L M . The squares \A p 



Nh', I = e, fi) are the light-to-heavy mixings, and sl in pip means the 
2 of the penguin diagrams A p ,, for the parameter choices (|3Tj) and (|32[) . 



are approximately proportional to (s£) 2 (s v £) 2 , as shown in Appendix of Ref. [la ] on the basis of unitarity of matrix 
U . The quantity Br(M — > e + e~) is dominated by the penguin diagrams, as seen earlier in this Section. Further, the 
other quantities considered later on in this work (cross sections for e~/i + — * e~e + and for e~/i + — * WW), are also 
approximately proportional to (s^) 2 (s^ E ) 2 . 1 While each of the two mixings separately is restricted by inequalities 



(s£) 2 < 0.005 , {s U l) 2 < 0.002 



(33) 



there is a strong restriction on the product of the two mixings: s^s^ 1 < 0.00012. The latter restriction is obtained 
in the two models, with parameter choice (|31[) and (|32p . when the models are required to give predictions for the 



For the box diagrams there can be substantial deviations from such proportionality, especially when matrix U is not real. 
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branching ratio Br(/i — > e 7) which respect the experimental upper bound [251 ]: 

Br(jT -> e~ 7 ) < 1.2 x 10 _u . (34) 
More specifically, the aforementioned upper bound for s v £ is obtained from restriction (]34p because 

Br( Ai -^e- 7 )=a em |-|G^| 2 , (35) 

and it can be shown that |G^| 2 < {l/4)(s£ 1 ) 2 (s l £) 2 , where approximate equality is obtained when the parameters 
are chosen according to Eqs. (f3"Tj) and (f3"2")l in models I and II, respectively. We will choose s v £ = s£ = s^. Therefore, 
we will take in (|3"Tj) and (f3"2")l the saturated values 

(sl) 2 = = i^f = 1-2 x 10- 4 . (36) 

For simplicity, we restricted Model II to cases with no CP violation (8\ = S 2 = 0). The heavy neutrino masses Ml and 
M2 are restricted to be above 100 GeV, and we take the convention Mi < M2. There is also an upper bound for Mi 
and M 2 in the form M 2 ^2 e \Bih\ 2 < ^M^/aw, with £ = e, /j, @, due to the breakdown of perturbative expansion to 
one loop above this bound. In Model I, this bound, in conjunction with restrictions (I36p . translates into a restriction 
which is approximately of the form Mi • M2 < (50 TeV) 2 , while for Model II the bound allows any value of M2 (> Mi) 
as long as Mi < 60 TeV. 

The numerical results for the branching ratio Br{M — > e + e~) defined in Eq. ([28]) are given in Fig. [3] and Fig.|H for 
Model I and Model II, respectively. We see that the branching ratios are extremely small, mainly because of the factor 
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FIG. 4: Same as in Fig. [3] but this time for Model II, for the (approximately) optimal choice (132 I t for the Dirac mass parameters. The 
curves, including the ones with M\ = 20 and 50 TeV, are in this model allowed by the PUB. 

~ 10~ 12 appearing in Eq. (|28p. Two major sources for this suppression are a factor cfi em coming from the muonium 
wavefunction and a factor a 2 m from the one-loop character of the squared amplitude. Another source of suppression 
are the small upper bounds (|36p on the mixing-parameters. From Figs. [3] and 0] we see that the branching ratios in 
in both models canot surpass 10 -19 . These values are so small that they seem out of the reach of experiments in the 
foreseeable future. 



/' 



+ e" -> e+e~ COLLISIONS AT HIGH ENERGY 



In addition to muonium decay, one can study the free scattering process /i + e — > e + e at higher energies. The 
cross section in the center-of-mass system (CMS) for this process at energies well above the muon mass follows from 
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the same amplitudes (fl2|) . (|T3|) and (fl~4)) . but the form factors enter in a quite different combination because of the 
different kinematic regime. Indeed, for CM energies y/s in the range m M -C y/s <C Mw, the cross section is: 



e e 
1 



12tt V4M^ 



(37) 



(1 - 2 sin 2 



)F e / + 2sm 2 6 w F^ + ^Z e \ 



4sin 4 ^|F7-F^| 



The numerical results, for Models I and II at various chosen heavy neutrino masses Mi and M2, are given in Figs. [5] 
(a) and (b), respectively. The elements of the Dirac mass mo are again chosen in the optimized forms (|31[) and 
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FIG. 5: The cross section for the LFV scattering process e~ fi+ — > e~ e + as a function of CMS energy ^/s in (a) Model I, and (b) Model 
II. In all cases, each pair of curves corresponds to a given value of Mi and to the values of M2/M1 = 1 (solid line) and M2/M1 = 2 (dotted 
line). 



respectively, and the mixing parameters have the upper bound values (|3l)|) . 

As shown in Figs. [3 the cross sections in both models are of the same order of magnitude, the ones in Model I being 
slightly larger than the corresponding ones in Model I. These cross sections are however very small. If one expects 
integrated luminosities not much larger than a few fb _1 at an hypothetical muon collider, then an experimental test 
of this process is not likely to be viable. 



VI. UNSUPPRESSED SCATTERING CHANNEL efi -> WW 

To complement our analysis of the LFV scattering process fi + e~ — > e + e~, we consider the case of open W produc- 
tion. 

When we consider e~/i + scattering at CMS energies above the W + W~ production threshold, the unsuppressed 
tree-level channel e~/i + — > W~W + opens up. This process is mediated by a neutrino exchange, as depicted in 
Fig. [BJ Here we will limit ourselves to estimate the cross section of the process. However, it is understood that the 
experimental observation of this process is hindered not just by the smallness of the cross section, but also by the 
rejection of background events that contain neutrinos in the final state. 

The amplitude T for fj, + e~ — > W + W~ is 

? = -^B^BIju rrr T--1 - h)KP2)*PLu{k 2 )] , (38) 

I |_(pi - kiY - nij + ie\ 

where e is the W-polarization vector. Consequently the total cross section for the process is: 

, 2 \ 2 N L +N H 

a(e-^^WW;s) = — -j ]T ( ir s 1^ H, , n;„) A„ (39) 

V J i,i=l 

= 64~ (10 h £ {BtjB^jBejBlj) (Ajj + A U A u A u ) , (40) 

W ^ ' S I,J>N L 
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where the explicit expressions of the coefficients Aij(mf, m?, Myy, s) are written in Appendix. To go from Eq. ([39]) 
to Eq. (|40p we used unitarity of the mixing matrix (the index t in Ay denotes a light neutrino state, whose mass we 
can neglect). The Dirac mass matrix mr> for models I and II is again chosen in the forms (|31|) and (|32|) . respectively, 
and the mixing-parameter values are those of Eq. (|36[) . 




FIG. 6: The tree-level process e /z+ — > W W + via exchange of a Majorana neutrino. 

The numerical results for Models I and II, and various chosen heavy Majorana neutrino masses Mi and M2, are 
given in Figs. 0(a) and (b), respectively. 




165 1000 10000 50000 165 1000 10000 50000 

s 1/2 [GeV] s 1/2 [GeV] 

FIG. 7: The cross section for the LFV scattering process e~ fi + — > W~W + as a function of CMS energy yfs in (a) Model I, and (b) 
Model II. Each pair of curves corresponds to the indicated value of Mi, and M2/M1 = 1 (solid line) or 2 (dotted line). 

We can see that the open W + W~ production cross section at y/s ~ 200 GeV is between 10 -5 and 10~ 4 fb, which 
is in general several orders of magnitude larger than the cross section for e + e~ production for CMS energies below 
the W + W~ threshold, 2 but it is still very small. As a reference, the point cross section a(e + e~ — > 7* — > at 
a/s = 200 GeV is close to 2 pb, which is 7 orders of magnitude larger. 

VII. CONCLUSIONS 

In the present work we have studied the possibility to observe lepton flavor violating effects in /i + e~ annihilation 
processes induced by heavy neutrino scenarios. We use two scenarios: one that is a plain seesaw model that violates 
both lepton flavor and lepton number (Model I) and another that has a richer spectrum of neutral fermions and violates 
lepton flavor while conserving lepton number (Model II). In a previous work[l6| we have studied these scenarios in 
the case of muonium to antimuonium conversion, which is box-diagram dominated and reflects the underlying process 
/i + e~ —>■ [i~e + . Here we were interested in the process /i + e~ — > e + e~ which is in general penguin-diagram dominated. 



2 If Mi and M2 are very large (> 10 TeV), the two cross sections at ^/s ~ 200 GeV are comparable, ~ 10 4 -10 5 fb. 
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We studied this process both at low and high center-of-momentum energies. At low (non-relativistic) energy, this 
corresponds to the muonium decay M — > e + e~. In both scenarios the branching ratios are exceedingly small (less 
than 10~ 19 ), basically because, besides the smallness imposed by a LFV process, there is an extra factor a^ em arising 
from the muonium wavefunction. To avoid this suppression, the alternative is to consider the collision /^ + e~ at higher 
energies. We find the cross section for /i + e~ — * e + e~ to rise with energy, but reaching at most 10 -5 fb at y/s ~ 50 
GeV. If one expects integrated luminosities not much larger than a few fb -1 at a future muon collider, then this process 
is not likely to be observed. Finally, we have pushed to very high energies, namely above the W + W~ production 
threshold. Only in that case we find that the two scenarios, Model I and Model II, would induce cross sections for 
/i + e~ — > W + W~ up to the order of 10 _1 fb or 1 fb. This process, however, has the additional experimental difficulty 
that it needs to be distinguished from the standard background where two neutrinos are present in the final state. 
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APPENDIX A: FUNCTIONS A tj FOR e/i -> WW SCATTERING 



The coefficient functions Ay (mf , to 2 , Myy> s) appearing in the sum (f39|) are long expressions obtained by calculating 
a Trace that follows from squaring the amplitude (|38|) . For this, we used the program Tracer [23]. When neglecting 
m M in comparison to Mw and y/s, the result for i ^ j is 



Ay 



[48 (mi - rrij)(mi + n%j)M^r 



|(to; — mj){mi + rn,j)sK 



Ajj — 



x [ - 6(m- + to 2 to 2 + m|) + 6(to 2 + m 2 )M^ - 24M^ - (3(m 2 + to 2 ) - 20Af^)s + s 2 
-6 [(to 2 - Ml/f(m\ + AM^) + (m? - 2to 1 M 2 / ) 2 s] A 
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where we used the short-hand notations 
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